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Abstract 

In this work, the delay limited capacity (DLC) of orthogonal frequency division multiplexing 
(OFDM) systems is investigated. The analysis is organized into two parts. In the first part, the impact of 
system parameters on the OFDM DLC is analyzed in a general setting. The main results are that under 
weak assumptions the maximum achievable single user DLC is almost independent of the distribution 
of the path attenuations in the low signal-to-noise (SNR) region but depends strongly on the delay 
spread. In the high SNR region the roles are exchanged. Here, the impact of delay spread is negligible 
while the impact of the distribution becomes dominant. The relevant asymptotic quantities are derived 
without employing simplifying assumptions on the OFDM correlation structure. Moreover, for both 
cases it is shown that the DLC is maximized if the total channel energy is uniformly spread, i.e. the 
power delay profile is uniform. It is worth pointing out that since universal bounds are obtained the 
results can also be used for other classes of parallel channels with block fading characteristic. The 
second part extends the setting to the broadcast channel and studies the corresponding OFDM DLC 
BC region. An algorithm for computing the OFDM BC DLC region is presented. To derive simple 
but smart resource allocation strategies, the principle of rate water-filling employing order statistics 
is introduced. This yields analytical lower bounds on the OFDM DLC region based on orthogonal 
frequency division multiple access (OFDMA) and ordinal channel state information (CSI). Finally, the 
schemes are compared to an algorithm using full CSI. 
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Index Terms 
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I. Introduction 

Multiple degrees of freedom in fading channels allow reliable communication in each fading 
state under a long term power constraint. This is due to the possibility of recovering the 
information from several independently faded copies of the transmitted signal. The rate achievable 
in each fading state is called zero outage capacity or alternatively delay limited capacity (DLC) 
[1]. Not only multiple input multiple output (MIMO) channels but also frequency selective 
multipath channels offer multiple degrees of freedom. This is in contrast to single antenna 
Rayleigh flat fading channels, where a DLC exists only if zero is not in the support of the 
fading distribution [2]. Since the DLC does not involve a decoding delay over multiple fading 
blocks if the variation of the fading process is slow enough, it can be considered as an appropriate 
limit for delay sensitive services, which become more and more important recently. Unlike the 
DLC, the traditional ergodic capacity strongly depends on the correlation structure of the fading 
process and generally implies an infinite decoding delay. 

This work investigates the DLC of frequency selective multipath channels in the context of 
an orthogonal frequency division multiplexing (OFDM) broadcast (BC) channel. OFDM can be 
considered as a special case of parallel fading channels with correlated fading process. Pioneering 
work on this topic was carried out in [3] where the general single user outage capacity was 
investigated of which the DLC is a special case. The optimal power control law is derived 
which is such that bad channels below some threshold are simply switched off. In [4] the limiting 
performance in the high signal-to-noise (SNR) regime of the DLC of parallel fading channels 
with multiple antennas was characterized assuming that the fading distribution is continuous. 
In [5] the impact of spatial correlations was studied. Further work was carried out in [6], [7]. 
Unfortunately, these results do not carry over to the OFDM case: since the subcarriers are highly 
correlated due to oversampling of the channel in the frequency domain the fading distribution 
is commonly degenerated which significantly complicates the analysis. This particularly affects 
the critical impact of the delay spread and the number of subcarriers. Hence, the behavior of 
the OFDM DLC is not clear yet. 
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Our main contributions are the following: First, we analyze the impact of system parameters 
such as delay spread, power delay profile (or the multipath intensity profile) and the fading 
distribution for the single user DLC in a general setting. We focus on two cases: the behaviour 
at high SNR and at low SNR. Such approach has been frequently used in the analysis of 
channel capacity even if the capacity itself is not completely known [8], [9]. The low SNR 
regime is characterized by its first and second order expansion. It is shown that to become 
first order optimal it is sufficient to serve only one of the best subcarriers regardless of the 
fading distribution. The corresponding limit is almost independent of the fading distribution. 
The second order limit is also calculated and depends generally on the number of supported 
subcarriers, i.e. when the fading distribution contains point masses. The quantities are shown 
to exist for a large class of fading distributions and are explicitly calculated in terms of the 
delay spread for the Rayleigh fading case. For the high SNR regime several universal bounds 
are calculated. These bounds culminate in a convergence theorem that generally characterizes 
the high SNR behaviour under very weak assumptions. Most important, there will be no need 
for concepts like almost sure convergence of the empirical distributions etc. (as used in [4]) and 
one approaches ergodic capacity relatively fast (the difference decreases with order log -1 (L) in 
a channel with L uniform taps and Rayleigh fading) even if the fading gains are not independent. 
The corresponding convergence processes are characterized. It is worth pointing out that these 
results not only hold for the OFDM case but also for other classes of parallel block fading 
channels. Finally, we provide a convergence result with respect to ergodic capacity. 

In the second part we focus on the broadcast scenario, complicating the analysis significantly. 
To begin, we present an algorithm which is capable of evaluating the OFDM BC DLC region up 
to any finite accuracy. This is a challenging problem, since for each fading state the minimum sum 
power supporting a set of rates has to be found [10]. To get a guideline for algorithm design, we 
subsequently derive lower bounds on the single user OFDM DLC based on rate water-filling and 
order statistics. These single user bounds are the point of origin for the development of simple 
analytical lower bounds on the OFDM DLC region based on orthogonal frequency division 
multiple access (OFDMA). The involved use of order statistics has a positive impact on the 
feedback protocol. In the low SNR regime, nearly the entire OFDM BC DLC region and in the 
high SNR regime a significant part of it can be achieved with these schemes without any form of 
time-sharing. Further, a practical OFDMA algorithm based on rate water-filling assuming perfect 
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CSI is introduced. This scheme outperforms the bounds based on partial CSI and might serve 
as a benchmark for other OFDMA minimum sum power algorithms. The results are illustrated 
by simulations. 

The remainder of this paper is organized as follows: Section II introduces the OFDM system 
model. In Section III and IV, the behavior at low and high SNR is studied in detail for the single 
user case. Section V contains the characterization and computation of the OFDM broadcast 
channel DLC region. Subsequently, in Section VI lower bounds on the OFDM BC DLC region 
are derived. We conclude with some final remarks final Section VII. 

A. Notation 

All terms will be arranged in boldface vectors (where m refers to users, k to subcarriers, I to 
path delays as a guideline) and the corresponding indices will be omitted if there is no ambiguity. 
Common vector norms (such as for the Zi-norm) will be employed. The expression z ~ 
CM (0, 1) means that the random variable z = x + jy is complex Gaussian distributed, i.e. the 
real and imaginary parts are independently Gaussian distributed with zero mean and variance 
1/2: x,y ~ A/"(0, 1/2). The expectation operator (e.g. with respect to the fading process) will 
be denoted as E (respectively E h ). Pr(A) denotes the probability of an event A. We write 
f (x) ~ g (x) if f(x)/g(x) — > 1 if x — > (or x — > oo) and all logarithms are to the base e 
unless explicitely defined in a different manner. 

II. System Model 

Assume an OFDM broadcast channel with M users from the set M. :— {1, ... , M} and K 
subcarriers from the set /C := {1, . . . , K}. The sampled frequency response of each user is by 
means of Fast Fourier Transform (FFT) given by 

L 

V— V _ 27Tj(i-l)(fc-l) 

h m ,k = e K , k e K (1) 

1=1 

where L < K is the delay spread and c m j are the complex path gains which are i.i.d. according 
to c my i ~ CM(0, (Jm,i)- The vector cr m = [t7 m ,i, cr m ,K] T is called the power delay profile (PDP) 
of the mth user's channel. The variances a my i are assumed to be strictly positive for all m and I 
and the channel energy is normalized Htr™^ = 1 for all users m G M. We say that the channel 
of user m has a uniform PDP if a^ m = . . . = ul,™ and a non-uniform PDP otherwise. Note that 
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in practice the PDP is typically non-uniform. Furthermore, the channel gains are not spread over 
the entire frequency band. Then, our results hold approximately and serve also as a performance 
limit. The channel (path) gains are defined as h myk '■= \h m ,k\ 2 (respectively c m ,i :— |c mi i| 2 ). The 
distribution of the channel gains is called the (joint) fading distribution. In case of complex 
Gaussian distributed path gains, i.e. c k ~ CM (0,1/ L) the channel gains follow a exponential 
distribution with Pr (h k > x) = e~ x . This case corresponds to Rayleigh fading. 

Let x my k with power p myk be the signal transmitted from the base station to user m on carrier 
k and let x m = {x mjl , x m , K } T ,Pm = \p m ,i, -,Pm,K] T the stacked vector of transmit signals 
and transmit powers for user m, respectively. Further assume that the system is limited by a sum 
power constraint E ^X)m=i ll*™!^) — -^-P* where P* is the power budget. Then the system 
equation on each subcarrier k can be written as 

Vm,k h rn<k ^ ^ %i, k ~l~ n m,ki k G /C, 

where y m>k is the signal received by user m on subcarrier k, and n m>k represents (without loss of 
generality) normalized noise n m>k ~ CM (0, 1). Let us define a decoding order n, such that user 
7r (M) is decoded first, followed by user n (M — 1) and so on. Assuming ideal superposition 
coding at the transmitter and successive interference cancellation with the decoding order it at 
the receivers, the rate of user 7r(m) is then given by 

= 1 E lQ g ( 1 + i h K(m) '% Am) ' k ) in [nats/Hz/s]. (2) 

K t^l V 1 + h < m )> k ^n<m PMn),k J 

For an instantaneous channel realization h = [hi : i, hi : K, h 2 ,i, 1im,k} t (respectively c = 
[ci,!, C\ : l, c 2; i, cm,l] T is the vector of path gains) the OFDM broadcast channel capacity 
region is given by the union over all power allocations fulfilling the sum power constraint P* 
and over the set of all possible permutations IT: 

C BC (h,P*)= |J [H:R Am) <R n(m) ,meM] 
Tren 

Here, R = [R\, Rm] T and R^m) is the rate of user n (m) defined in (2). 
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III. OFDM SINGLE USER DLC 

A. An implicit formulation of the single user OFDM DLC 

First let us define the single user DLC Cd for OFDM, which is a special case of parallel 
fading channels. The user index m is omitted in this section. 

Definition 1: A rate R* is achievable with limited delay (zero outage) under a long term power 
constraint P* if and only if for any fading state h there exists a power allocation p' solving 

min jlp'lli subj. to R* < R(h, p') (3) 

and 

Eh(Hp'Hi) <P*. 

Furthermore R* is called the delay limited capacity C d (P*) if and only if 

Eh(|| P '||i) = P*. 

Obviously, the delay-limited capacity is smaller than the ergodic capacity under the same 
average power constraint because the requested rate has to be supported in every fading state. 
Observe that we choose to state the DLC region as a definition rather than a theorem as done 
in [1]. There it is shown that our definition coincides with rates that guarantee arbitrary small 
erroneous decoding probability (dependent on the coding delay) for any jointly stationary, ergodic 
fading process when the codewords can be chosen as a function of the realization of the fading 
process. This can be interpreted as a performance limit for traffic where the information has to 
be delivered within one coding frame for almost all fading realizations and the channel varies 
slowly enough such that the transmitter can track the channel. Moreover, by our definition we 
see that since the sum power minimization problem has to be solved in every fading step the 
DLC represents also a performance limit for practical minimum sum power algorithms. We will 
present an example in Sec. VI. 

The single user DLC was already examined by several authors in the context of systems with 
multiple antennas and parallel independent fading channels [4], [5]. Nevertheless, we re-derive 
the DLC here using the principle of rate water-filling, which leads to an interesting perspective. 
This characterization of the DLC will be used later on to derive lower bounds on the OFDM 
broadcast channel DLC region in Section VI. Denoting the rate on subcarrier k as r k , the 
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optimization problem for each fading state in (3) is equivalent to 

K K K 

min ^ Vk ~ A ^ r k subj. to ^ r k > KC d . 

k=l k=l k=l 

Using the relation between power and rate on each subcarrier this can be expressed as 

min^ ( ( —— — - - \r k 



k=i 



where A > is to be chosen so that J2 k=1 r k = KCd- The resulting optimality conditions are 
given by 

= Vfc, J2 rk = KCd {A) 



hk 



k=i 



with [•]" = min{-,0}. Note that by taking logarithms and solving for r k eqn. (4) can be 
interpreted as a water- filling solution with respect to the rates r k . Combining the optimality 
conditions for all K subcarriers and solving for the Lagrangian multiplier A yields 

exp ( °f-) 

llfcez>(h) ^fe h 

where the random variable V (h) C K. denotes the set of active subcarriers and d h := \V (h)|. 
Note that the allocated power is given by p k = A — h^ 1 for any k e V (h) and zero otherwise. 
Substituting (5) in the average power expression given by 

(6) 



s k=l 




with [-] + = max {•,()} we obtain the single user OFDM delay limited capacity d with power 
constraint P*\ 



Since the denominator in (7) can be bounded by a constant, the delay limited capacity C d is 
greater than zero if and only if 

1 —^ T dF h (h) < oo. (8) 

llfce©(h) h k 

Here, F h denotes the joint fading distribution function. The class of fading distributions for which 
(8) holds is called regular in [4]. It will become apparent in the following that the correlation 
structure of the channel gains in OFDM provides the main challenge in proving and analyzing 
regularity according to (8). 
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B. Suboptimal power allocation strategies 

Let us introduce a suboptimal power allocation for the single-user case that is used in VI. 
It is evident from the expression for the single user DLC that the major difficulty is the rate 
water-filling operation for each channel realization. To circumvent this difficulty which results 
in a prohibitive complexity for the multi-user case we introduce the notion of rate water-filling 
for average channel realizations. The idea is to use simply the information which subcarrier is 
the best, the second best and so on and to allocate fixed rate budgets to the in that way ordered 
subcarriers. For the analysis we need the following definitions. For a given vector h of real 
elements let us introduce the total ordering 

hk[K] > h k [ K -i] > > h k [i], 

i.e. h k [i] is the minimum value and h k [M] is the maximum value. If h is a random variable then 
the distribution of h k [ p ] is known to be the p-th order statistics. The the p-th order statistics can 
be explicitely given for K independent random variables h k with distribution F and density / 
from standard books. The p-th order density is given by 

f hm (x) = Kf (x) ^ F^ (x) (x)) K ~ p (9) 

Based on the order information and the distribution we can now deduce a fixed rate allocation on 
the subcarriers, avoiding the water-filling procedure in each fading state. The idea is to allocate 
a fixed rate budget to the p-th ordered subcarrier. Defining the terms 

poo 

C P = / h' 1 dF hk[p] (h) 
Jo 

and using these factors the optimal rate allocation is now given by solving the optimization 
problem 

e R k\v] - l 

mm — — XR k , yR k ip]=KC d , 

R kip]>° % |— ; 

which we have already solved by rate water-filling given by 

K 

= 0, ^2 Rk{ P ] = KC d - 
p=l 

Obviously, this suboptimal scheme has some impact on the feedback protocol: Since only the 
order statistics are exploited, it seems to suffice to feed back the ordering of the subcarriers. This 
affords much less feedback capacity than perfect channel knowledge would require. However, 



log (e) 



3 -Kfc[p] 



Cr 1 



- A 



8 



it is not straightforward to translate rate water-filling to a power allocation strategy, since for 
achieving a certain rate the channel has to be perfectly known. Allocating fixed power budgets 
is possible. The consequence is, that since the p-th order channel is a random variable, mutual 
information becomes a random variable once again not guaranteeing a certain rate in each state. 
However, the variance becomes much smaller. 

There is an interesting second power allocation where the powers asserted to the subcarriers 
are all the same. It is easy to see that the allocation according to 

Pk = — i— K tj^tt, keIC, (10) 

E h (nf =1 h~ k 1/K ) 

always leads to a rate higher than the requested rate with equality in the high SNR region. 
Hence, this is also a suboptimal solution. The bounds are illustrated in Fig. 1. 

It is of great interest to understand the impact of the delay spread L and the power delay 
profile er as well as the fading distribution on the OFDM delay limited capacity. In case of the 
OFDM broadcast channel, an analytical characterization is nearly impossible. Thus we carry out 
an analysis for the single user OFDM DLC in the following. Note that this matches the behavior 
on the axes of the OFDM BC DLC region where only one user is active. Hence the results give 
insights for the broadcast case as well. Since the expression in (7) is still very complicated, we 
focus on the behavior in the low and the high SNR regime and carry out a detailed analysis. 

IV. Impact of system parameters 

A. Scaling in low SNR 

1) Impact of delay spread and fading distribution: First, we characterize the first and second 
order behavior at low SNR. For ease of notation we define := llhll . 

II 1 1 oo 

Proposition 1: Suppose that E h (h^) < oo. Then, the following limit holds: 

o P* E h (/i-!) 

Proof: Our starting point is (7) where we use the McLaurin-expansion of the exponential 
function exp (x) — 1 + x + o (x) up to the linear term to obtain a lower bound on the required 
power P*. 

Fix now e > and use the following strategy: set the 'virtual' channel gain of any subcarrier 
k with hk > hoc — e to hoc. Denote by Xh (e) the multiplicity of the number of subchannels that 
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are assigned the maximum channel gain by this strategy for any channel realization. Using (7) 
we have for sufficiently small P* an upper bound on C d which is given by 

= C d E h f-L) (13) 



and we obtain 



Q £ w (14) 



for any e > 0. 

For the upper bound on P* fix e' > 0. This time set the number of supported subcarriers to 
one and support one of the set with maximum channel gain. Then we have for sufficiently small 

Hence we have 

C * £ (1 + /E,, (0 <16) 
for any e' > 0. Combining both the lower and upper bound yields the desired result. ■ 

This quantity also reveals the minimum energy per bit, at which reliable communication is 
possible under a limited delay. This is in analogy to [8], where this quantity was derived for the 
ergodic capacity of a Gaussian channel. Moreover, the lemma states that albeit it is generally 
suboptimal, serving one of the best subcarriers becomes optimal in the low SNR region. This 
can be easily seen since the multiplicity Xh(e) of the attained maximum channel gain vanishes 
in the expressions for the lower and upper bound. 

To gain more insights, the remaining sub-linear term defined by 

A d (P*):=C' d (0)P*-C d (P*) (17) 

is calculated next. The following proposition tells us that while for the linear term it did not 
matter if the distribution contains point masses it does matter for the sub-linear term: 
Proposition 2: The following limit for the sub-linear term holds: 

Here, Xh is the (random) multiplicity of subchannels with maximum channel gain. 
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Proof: Our starting point is again (7) where we now use the McLaurin-expansion of the 
exponential function up to the quadratic term, i.e. exp (x) — 1 + x + 0.5x 2 + o (x 2 ). 

By the same strategy as above fix e > and set the "virtual" channel gain of any subcarrier 
k with h k > hoc — e to h^. Then, we obtain for sufficiently small P* 

p* > Eh ( Cj+ y e ») . 

The equation is an upward open parabola in C d where one zero is negative and one is positive 
where the latter is increasing in P*. Solving this equation for C d yields the inequality 

Eh fc 1 ) 



C d (P*)< 



KE h (x- 1 (e) h-J) 



K(h-J) | 2P* 



(20) 



K 2 K (Xh 1 (e) h-j) KE h (x h X (e) h~j) ' 
Expanding the square root function yields 

Cd(p) - ^T) p ~w^mm {p } ■ <21) 

Subtracting the first order expression (11) from (21) we arrive for some e' > at 

A 2 (1 + 6) E h (fc-i) ^ 

" 2(l + e) E h (/,-!) (22) 
and thus have established a lower bound on A d (P*) for any e, e' > 0. The last inequality (22) 
follows from the following argument (which is frequently used in the sequel): observe that 
Xh (e) > 1 and, almost surely with respect to the fading distribution, for any realization h we 
have 

XhM^^Xh 1 ^ 1 , e^O, (23) 
and provided that E h (Xh 1 ^^, 1 ) < E h (h^-) < oo we obtain by dominated convergence [11]: 

limEh (xh 1 (e) C) = E h fa 1 /*" 1 ) (24) 
In analogy to the derivation of the first order behavior we get 

MP)£ — ^(o — (p) ' <25) 

Combining (22) and (25) leads to the desired result. ■ 
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Corollary 1: Since Xii( e ) > 1 a simple upper bound on the second order term in (18) is given 

by 



A, (P*) AT 

-p^o 1 (P*) 2 ~ 2 [E h (/^ 
Note that the bound from Corollary 1 is consistent with the result in [12] where it is shown 



limsup^^; — — < j- (26) 



that 

C,(p.)~llog(l + i ^), (27) 

This can be easily checked by differentiating the expression in (27) twice. 

Proposition 3: Suppose that the joint fading distribution is absolute continuous. Then, the 
following limit holds: 

, im Mn = k 

f-o ( P .f |E h (ft-i)] 2 
Proof: We have to prove that the set of events where the maximum is taken on by more 

than one subcarrier has probability zero. 

By the absolute continuity of the joint fading distribution (which is preserved under unitary 

mappings) this is equivalent to show that the set that contains all events where the maximum is 

not unique has Lebesgue measure zero. To see this we consider sets of the form 

{h G : [J {h n (i) = . . . = h n (i), h n ( i+ i), . . . , h^ K )} \ 
l<i<K ) 

where it is the permutation that yields an decreasing order among the channel gains 

Observe that any set has measure zero. Since there can only be 2 K such possible sets the union 
of these sets has measure zero as well. Now, as our regarded set is in the union of the constructed 
set, the set has also measure zero. ■ 

Hence, appealing to Prop. 1, 2 and 3 the forthcoming analysis reduces to the study of the 
expected maximum of the channel gains. However, the expressions do not show how the DLC 
depends on the system parameters which we investigate by means of an asymptotic analysis, i.e. 
for large L, K. This analysis turns out to be quite accurate even for very small L. 

Remark 1: It is important to note that for the asymptotic analysis we will let go L and K, 
K > L, to infinity which is indicated by the index n, i.e. K n , L n — > oo as n — > oo. We assume 
that for all L n the complex path gain vectors c n = [ci, . . . , clJ t are defined on the same 
probability space, i.e. to each c n G C Ln there is by means of (1) h n = [hi, . . . , Iik} t G C Kn . 
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While the distribution does not change for the L n _ x first random variables of the vector c n the 
distribution of the K n - X first random variables in h n might change since they depend on the 
FFT structure. We have to keep this in mind for the forthcoming analysis. 

It was shown in [13] that the (continuous) maximum value of the frequency response equals 
log (L) with large probability even for moderate L for the following distributions: 

Cj independent and iid (both in real and imaginary parts) Vi, 

E h (Ke 2 (50) = .5/L,E h ( e ^ e ^) = e -^-^ + El 3 a^+o(^) ^ (29) 

for all \u\ < d, d > 0, a 3 , a 4 , a 5 G C 

Here, (c) denotes the real part of the complex number c (Sm is the imaginary part). Note that 
the condition on the characteristic function of the real part of the path gains implies finiteness 
of the moments up to order six of the corresponding distribution. The following theorem proves 
that the log (L) result holds also for the maximum of the sampled frequency response regardless 
of the sampling set. Furthermore, it reveals the exceptional role of uniform PDP. 
Theorem 1: Suppose that the fading distribution belongs to JF L then we have 

Pr (log (L) — g (L) < hoc < log (L) + g (L)) = 1-0 (log" 4 (L)) 

with g (L) = 4 log [log (L)] for large L and arbitrary K > L. Furthermore, the upper bound 

Pr (/loo > log (L) +g(L))<0 (log" 4 (L)) 

also holds when the PDP is non-uniform. 

Proof: see Appendix A. ■ 
We can apply this result to the DLC assuming uniform PDP where we have to show that 
from the convergence in probability given in Theorem 1 follows convergence of the expected 
maximum of the channel gains. This can be achieved if the set of distributions is somewhat 
more restricted compared to (29) in the sense that the their behavior in the neighborhood of 
the zero is "sufficiently well". By this we mean, that the distribution function of c\ is Lipschitz 
continuous in an e-neighborhood of the zero. Let us denote this class of fading distributions by: 

^Lo (k s ) ■= {the conditions (29) on q hold Vi 

h e is generated from c e by means of (1) 

F C1 (x) < k s x, < x < e, k s , e G 
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The Lipschitz continuity is essential in the next Lemma. 

Lemma 1: Suppose that the fading distribution belongs to jF Lo (k s ) then the following limit 
holds for sufficiently large L and arbitrary K > L: 

Jim Q ^P = log (L) + O (log [log (L)]) 

Furthermore, the DLC is maximized (with respect to the leading order term) by uniform PDP 
("order-optimal"). 

Proof: The last statement follows from Theorem 1 and, hence, we assume uniform PDP. 
By the same theorem there exist constants 7, k e IR ++ so that 

Pr (log (L) - 7 log [log (L)\ < < log (L) + 7 log [log (L)\) 

K 



> 1 - 



log 7 (L) 

for sufficiently large L. Setting e_ = 1 — e and e + = 1 + e where e := 7 log [log (L)] / log (L) 
the expectation can be written as 

Eh fc 1 ) 

= E h (/I" 1 ] h,, e [e_ log (L) , e+ log (L)}) Pr (h^ G [e_ log (L) , e + log (L)}) 

+ E h (/I" 1 ] ^ < e log (L)) Pr (^ < e log (L)) 

+ E h ( h-J \h OQ >e + log (L)) Pr (h^ > e + log (L) ) . 

Note that the crucial part is to derive an upper bound on the conditional expectation in the 
second term on the RHS of the last equation, i.e. when the maximum is small. We will show 
now that the conditional expectation is bounded too. Proceeding with the standard inequality 
E {\X\) < 1 + Y^=\ P r (1-^1 > n ) an d using the "trick" that we can apply it to the conditional 
probability measure as well, the second term is given by 

Eh ( C| ». < «- log (L)) < 1 + £ Vl[{h ^t^^ m ' 

Using the inequality Hc^ = Hh^/if < h^, and < Hc^, define P (x) := Pr ({hoc < x}) 

and we obtain by independence of the path gains 

E h (/i- 1 |/i 00 <e_ log (L)) 

i+LM p/i\ 1 +00 /l , . \ 

e ?dik + 5i^ii 1 ft (n{^7}j (30 > 
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where [k s \ is the greatest natural number below k s . Hence 
RHS of (30) = 2 + [A; S J 



P(c_ log (L)) 



< 2 + [A; S J + 



k 



i=2+[k s ] 



£ * (t) 

Oi I t, I \ / 



E 7l - 



Pfe.logfL))^.,. 



(31) 



In the second step we assumed Lipschitz continuity of the path gain distribution function F Cl 
with Lipschitz constant below or equal k s , i.e. F C1 (x) < k s x. The series in (31) can be upper 
bounded as follows: 

1 



i=2+[k s \ 



E i< E 

i=2+[k : 
r+oc ^ 

A+ifc,i 



dx 



x 1 



dx 



< 



i+LM 

-L + l L Jl +LM 
1 



L - 1 



fc: 



-L+l 



1 



Hence, we obtain finally 

RHSof(31)<2 +W + p(£log(i))i _ 1 . 
The last term is finite for L > 1. Hence, we have the upper bound 

Eh (h'J) < ej log- (L) + (2 +i fc ;J) « + _A_ + c; i jog-i ( L) K 



log- (L) L — 1 to v 7 log- (L) 

_ f -i, m -i/ n A | (2+lfc 8 ])e_K A; s e_log(L) e_K 
" e - l0g (L) I 1 + log-- 1 (X) Z - 1 + 

= el 1 log" 1 (L) f 1 + O 



e + log- (L) 



1 

v7~ 



log-- 1 (L) 



for 7 > 1. Therefore, we have finally 



C d (P*)>-\og\l + KP\_ log(L) 



= -log(l + KP*e_log(L) 



1 + 
1-0 



log 7 " 1 {L))\ 
1 

log 7 " 1 (Z) 



A lower bound is obviously 



E„ (/»-) >e;Mog- (i ) (1-1^%)) 
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and we have 

C d (P*) < j£ log ^1 + KP*e + log (L) ^1 



1 log! l + KP*e + \og(L) 



log 7 (L) 
I -0 



log^ 1 (L) 



and the result follows. ■ 
We can conclude from the proof as follows: 
Corollary 2: The DLC is finite if there are either 

• at least two independent channel gains in the frequency domain or 

• at least two independent path gains in the time domain 

with Lipschitz continuous marginal distribution function in the neighborhood of the zero. 

Interestingly, the DLC compares favorably by the factor log (L) with the capacity of AWGN 
in the low SNR regime. Hence, the delay spread governs the DLC in this region. 

Let us make the bound explicit for the Rayleigh fading case. Note that due to the sum 
of independent complex path gains the samples of the frequency response are approximately 
complex Gaussian distributed anyway (however, the exact distribution would be very difficult to 
evaluate) [14]. 

Lemma 2: Suppose L divides K and under the assumption of complex Gaussian distributed 
path gains with uniform PDP the maximum channel gain is enclosed by the inequalities 

Pr (log (L) - 7 log [log (L)] < < log (L) + 7 log [log (L)]) 
>1 



log 7 (L) 

where 7 > and k > L _ lo ^—,^ and 7, L sufficiently large such that the bound makes sense. 
Proof: see Appendix B. ■ 
We can apply the result again to the DLC. 

Corollary 3: Under the assumption of complex Gaussian distributed path gains with uniform 
PDP the low SNR DLC is enclosed by 

1 log (1 + k,KP* log (L)) < C d (P*) < 1 log (1 + k 2 KP* log (L)) 

where 

7 log[log(L)]\ / I.Ik 



«i := 1 



log(L) A log 7 ^) 



16 



and 

/ 7 log [log tm i ijk \ 

for any 7 > 1, k > 0, and L > 1 from Lemma 2. 

Proof: The result follows if the constants in Lemma 2 are used in the proof of Lemma 1 
and by simple numerical check. ■ 
Remark 2: If L does not divide K it seems very difficult to get results that are asymptotically 
tight. However, good bounds are easily obtained by observing that the frequency response cannot 
arbitrarily overshoot between the samples. In fact, for the upper bound we can use a grid K' = 
aL, say by a factor a > 1, and by collecting these samples in h( a ) we have [15] 

llhll < ||h<°>|| cos-^—V 
11 Hoc - || Moo \2a) 

For the lower bound we set a = 1 and have 

||h|| > llhWII cosf— V 

II Moo - II Hoc \2K J 

Since we have just derived bounds for [[h^U^ , a > 1, we can tackle also the general case. 

An illustration is shown in Fig. 2 where we calculate (27) for different but low SNR. It is 
observed that the approximations are quite accurate for small L. The behavior of the OFDM 
DLC at low SNR and the corresponding first and second order approximations are depicted in 
Figures 4 and 5. It can be seen that the region where the approximations hold diminishes as 
the number of degrees of freedom increases. The bounds can be used to roughly estimate the 
performance of e.g. a cellular system at the cell border. Even though we have made not effort to 
optimize the bound we found by simulations that for Rayleigh fading with small delay spread 
the error is within reasonable span of the optimal curves. However, it is worth noting that for 
large ratios of L and K the range where the bound makes sense becomes increasingly small. 

2) Impact of power delay profile: The impact of the PDP has been touched already in Lemma 
1 proving the "order-optimality" of uniform PDR Let us now investigate the general case. The 
following expression can be used to get a bound for arbitrary PDR 

Proposition 4: For sufficiently low SNR an upper bound on the DLC is given by: 

1 / if(i+<*) p* \ 

c ' (n -^( 1+ i(wr)j (32) 
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Here, a £ (0, 1) is a global parameter that can be numerically optimized. The lower bound is 
independent of the order of the elements of the PDP and concave, thus Schur-concave. 

Proof: The proof follows from Prop. (4) and the inequality chain Hhj^ > || h|| ^ > 
K" 1 Uhllj = \\c\\ v Since both lower and upper bound are tight only in identifiable special 
cases there is some a £ (0, 1) that can be numerically found. The Schur-concavity is obtained 
from Prop. B.2. in [16, pp. 287] since Hcjl^ 1 is symmetric and convex and the expectation 
E(||c|| ] ^ 1 ) is independent of the ordering of the PDP. ■ 
The Schur-concavity implies that for fixed L and normalized PDP the bound is larger if the 
elements of the PDP vector are more "spread out". If L is large the bound approaches the low 
SNR AWGN capacity time a factor K a for uniform PDP due to the strong law of large number, 
i.e. 1 1 c 1 1 -l — > 1 almost surely, and hence a is of order log [log (L)] / log (L). 

B. Scaling in high SNR 

Defining the quantity h := Ylk=i ^ was P rove d in [4] that by using the suboptimal 

power control law (10) C d (P*) > log (P*/E h (h)) provided that E h (h) < oo, i.e. for regular 
fading distributions. We can extend this result to an upper bound without using any simplifying 
assumptions on the fading distribution. 

Proposition 5: Suppose that E h (h) < oo. For sufficiently large P* the DLC is upper bounded 



that for any fading state h we set the values that are below e to e. In other words we do not 
allow "virtual" channel gains below e. 

Define h% := max{/i fc ,e}. Then, using (7) we have for sufficiently large P* 



Obviously, the second term grows without bound as e — > for many fading distributions such 
as Rayleigh fading. Furthermore the growth depends on P*. Let us bound this term as follows: 



by 




Proof: We can use the following strategy for an upper bound Cd- fix e > and suppose 




(33) 
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we have 



Eh (£>=f 4™ 



Clearly, the term e is related to P*. Since the maximum channel gain is at least e and the 
underlying optimal power control law is water-filling the above equation (33) is certainly true if 

e 

which is a rough estimation. Hence, we obtain 

/ 1 \ P* 

and finally for any e > 

P* (1 + i) 



C d < log 

Now observe that Ylk=i i^l) ^ ^ h an< ^ 



E h (nf =1 (k 



k> k 



fc=i 

Hence, by dominated convergence 
provided that 



vfc=i 



E h (/i) < cx). 

■ 

The proposition states that as long as E h (/*) < oo the DLC lies in some target corridor 
determined by E h (hj . The following proposition was proved in [4] where it is shown that, under 
appropriate circumstances, serving all subcarriers equally is sufficient to achieve the limiting 
performance. 

Proposition 6: Suppose that E h (hj < oo. If the joint distribution is continuous, then 



C d (P*)~ log 



p* 
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Hence, appealing to Prop. 5, and 6 it suffices to evaluate the term E h (K) in the high SNR 
regime. The problem of whether or not the high SNR quantity is non-zero is not touched upon 
in [4]. Let us therefore derive general conditions under which this is true. As before we could 
assume that the fading distribution is Lipschitz continuous on which is still too restricting 
though. 

Remark 3: Note that one is attempted to derive these conditions from the easier low SNR 
quantity since the DLC increases with SNR. However, this is in general misleading since the 
joint fading distribution might e.g. contain point masses on the boundary of the positive orthant 
rendering the high SNR term infinite while the low SNR term still provides a proper lower 
bound. 

The following proposition puts a general lower bound on the DLC and will be used in Lemma 
(3). In order to not overload the formalism we assume that the joint distribution on a equidistant 
subset with some distance a E N of b E N subcarriers possesses a density as follows. Define 

h (a, k x ) := [h kl+0a , h kl+a h kl+ab ] T , k x E (1, . . . , a). We define the following class: 

Fm(c s ) := {F h(aM) has density f h(aM y, 

/h(a,fci) an d all marginal densities are bounded by c s 
within 1/b- neighborhood for all k± E (1, . . . , a) 
h E is generated from c E by means of (1)} 

Observe that we have not excluded point masses in this definition. 

Proposition 7: Suppose that the joint fading distribution on this subset belongs to jF Hi (c s ). 
Then, for sufficiently large P* the following lower bound holds: 

C d (P*) > log (P*) - log (^bc s (jZT[) + t 1 " c *) 6 ) (34) 

If a can be chosen to be K/2, i.e. b = 2 for even K we have: 

C d (P*) > log (P*) - log (6c s + 2) (35) 
Proof: see Appendix C. ■ 
Clearly, the bound can be improved in case of independent subcarriers which will be used in 
Cor. (2). 
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Proposition 8: Suppose that the joint fading distribution can be written as Ph(a,fci) = lL=i / 
for all fci 1 , i.e. subsets of the subcarriers are independent, and that the marginal channel gain 
density / is finite everywhere. Then, for sufficiently large P* the lower bounds (34), (35) can be 
improved to give: 

b 



C d (P*)> log (P*)-Mog 
-log 



b-1 



\ V6 

£/(«r)-/(&r) -E^(/W)-/W)) 

i>l / i>l 



Here, < ^ /+ < &; /+ < oo are interval boundaries of supp (/) such that /' (h) < 0, h G 
[a", b~] and /' (/i) > 0, h G [a+, 6+] . If b is large and /' (h) <0,heR, with / (0) = 1 (e.g. 
Rayleigh fading), then (b/b - if -> e 1 and hence C d (P*) > log (P*) - 1. 

Proof: see Appendix D. ■ 
For the next proposition we need explicit the properties of the FFT. 

Proposition 9: Suppose that L, K are even and that the densities of real and imaginary parts 
of the complex path gain distribution fulfill /sR e (c i ),Sm(c l ) (x) < y^e - "^' 2 . Then, the following 
lower bound holds: 

a L TfL 1/2 / T \ L-1 N 



L L 



a 



C d (P*)>log(P*)-log 
Proof: see Appendix E. 
The latter proposition is universal but tailored to the Rayleigh fading case. More sophisticated 
bounding techniques can be obtained from mixing the techniques of Prop. (7) to Prop. (9) as 
discussed in the remark of Appendix E. In the following we assume without loss of generality 
that the marginal fading distributions are such that E h (log (h k )) is independent of k e /C (and 
finite) where K is supposed here to be even. 

Lemma 3: Suppose that E h (h) < oo. For sufficiently large P* the DLC is upper bounded by 

C d (P*)< log (P*) + H(F h ) + ^ 



'it is straightforward to see (by the structure of the FFT) that if the fading distribution is generated by a complex path gain 
distribution of which the density can be written as fg i (c») , c» £ C, where fg i is invariant under complex rotations, then the 
fading distribution is also invariant regarding k\. 
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where 



poo 

H(F h ) := / \og(h) dF hl (h) 
Jo 



and F hl is the marginal fading distribution. Furthermore, suppose that there is a sequence of 
fading distributions in jF Hi (c s ) such that K^ 1 J2k=i 1°S (^fc) — * ^ (-^h) in probability. Then, the 
bound is asymptotically tight, i.e. 

C d (P*)- log (P*) + H(F h ). 

Proof: Setting 

E h (h) = E h (exp [log (h)]) 

we get by Jensen's inequality 

E h (exp [log (h)])> exp (E [log (h)]) 

= exp ^ log (h) dF hl (hfj 

which is already the desired upper bound provided that H (F h ) < oo or equivalently E h (h) < oo. 

In order to show the tightness of the upper bound define the following random variable (i.e. 
partial sums): 

1 K 

^:=--5>g(/i fe )- (36) 

k=l 

Suppose that — > H (F h ) in probability. We have to show that 

£ h (exp 

-»■ exp (- log (h) dF hl (h)^j , K -> oo, 

which would follow if is uniformly bounded in K but is not true for the situation at hand. 
Using the set function !{■} and writing for some C > 

E h (exp (ftW)) = E h (exp (ft^) I {/j^ < C}) 
+ E h (exp p)l{/iW>C}) 
< E h (min {exp (h {K) ) , exp (C)}) 
+ E h (exp (/ l W)l{/ l W>C}) 
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K — > oo, 



yields for the first term on the RHS 

E h (min {exp (h {K) ) ,exp (C)}) 

-> exp / log (h) dF hl (h) , 
Jo 

since hS K ^ — > H (F h ) in probability and min {exp (/i^) , C} is uniformly bounded (and C 
sufficiently large!). Hence, we have 

limsupE h (exp (h (K) )) 

< H (F h ) + limsupEh (exp (/i (x) ) l{h^ > C}) 

n^oo 

= #(F h ) +limsup f hl{hS K) > C) dF h (h) 
Fix e > we have by the inequality shown in (48, Appendix C) 

K/2 



where VL C 
argument 



/ hl{h™ > C) dF h (h) < f[ / h\h} +K/2 dF h (h) 
{he : > C}. Then by geometric mean inequality and a "sandwich" 



hl{h (K) > C) dF h (h) 



h ! h f+K/2 ~ ( k f) 2 (K+K/2)) dF h (h) 



K/2 

<-Y 

i=i 
9 K / 2 

^e^Pr (hW>c) + C e 

where Pr (h^" 1 > C) — > for if — > oo (C again sufficiently large). The remaining constant 
C e > can be made arbitrarily small independent of K since F h e JF Hi (c s ) uniformly. On the 
other hand, since 

liminf E h (min {exp (h {K) ) , exp (C)}) > H (F h ) 

we have the desired result. ■ 
The required convergence in probability follows if: 
• either, the subcarriers are independent (or a subset), 
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• or the logarithmic channel gains are uncorrelated (or a subset) with 

1 K 

— ]TE h (log 2 (^)) -0, K^oo. (37) 

fc=i 

Interestingly, since H (F h ) < there is always a loss in capacity compared to AWGN. 
Furthermore, observe that the second statement (37) is substantially weaker than independence. 
It suggests that ergodic capacity can be achieved even if the subcarriers are not independent 
which is discussed in the next subsection. We can apply this result to the Rayleigh fading case. 

Theorem 2: Under the assumption of complex Gaussian distributed path gains, the upper 
bound 

C d (P*)< log (p*) + tf(F h ) + -i 

holds where 

PCX) 

H (F h ) := / log (h) exp (h) dh « -0.58. 
Jo 

The bound is asymptotically tight for uniform PDP and sequences (L n , K n ) where L n divides 
K n , i.e. 

C d (P*) ^ log (P*) + H (F h ) , n - oo. 

The convergence speed is given by: 

C d (P*) > log (P*) + H (F h ) + O (log- 1 (K)) 

(see eqn. (38) for constants) 

Proof: The first part of the theorem follows immediately from Lemma 3 and the fact that 
subsets of L subcarriers are independent. It remains to provide an explicit expression for the 
convergence speed. Let be defined as as in (36). Then by Lemma 3 we have to investigate 
the following terms 



E h (exp (fc<*))) 

< E h (min {exp (h {K) ) , C}) + E h (exp (h (K) ) I {h iK) > C}) 
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where C > H (F h ). Defining now the event A = { \h^ K) - exp (H (F h ))| < e} and its 
complement „4. c we have 

E h (exp (/j^)) < E h (min {exp (ftW) , C} I {A}) 
+ E h (min {exp (fcW) , C} I {^ c }) 
+ E h (exp (/ 1 W)I{^)>C}) 

The first two terms can be bounded as follows: since 

E h (min {exp (h {K) ) , C) I {A}) < exp (H (F h )) + e 

and 

E h (min {exp (/i (x) ) , C} I {A c }) <CPi (A c ) 

we need a bound on the probability Pr (A c ) dependent on e. The probability can be easily 
upper bounded by Tschebyscheff's inequality, i.e. 

Pr(A c )< " 



Ke 2 
where 

( x 2 =E h ([log(/i 1 )-E h (log^))] 2 ) 

and by choosing some sufficiently slowly converging zero sequence, e.g. tx = O (1/ log (K)). 
The third term can be upper bounded by observing that: 

E h (exp (hW)l{hW >C }) 

hl{h {K) > C} dF h (h) 

< ^ (h) 2 dF h (h) j 1 ^ (I {ftW > C}) 2 rfF h (h) 

< v / ^pO)(Pr({^ ) >C}))^ 

In the last inequality we employed Prop. 8 for the first integral. The probability can again be 
tackled with Tschebyscheff's inequality. It follows 



Pr ({/i<*> > C}) < 



a 2 



K(C-H(F h )Y 
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and hence: 

C d (P*) 



> log (P*) - (exp (H ( F *)) + ' K + ^- + — 



exp 1 / 2 (1) a 



/ 2 (C-H(F h )) 



An illustration is shown in Figure 3. 

1 ) Impact of power delay profile: The impact of the PDP has been touched already in Cor. 
2 showing its asymptotic optimality. Similar to the low SNR regime we have an upper bound 
describing generally the impact of PDP. 

Proposition 10: For sufficiently large P* (and K) the following upper bound holds: 

C d (P*)< log' 1 



EhdlcH- 1 ) 

The bound is independent of the order of the elements of the PDP and concave (thus Schur- 
concave). 

Proof: The proof follows immediately from the geometric mean inequality, i.e. for any h 
we have 

- K 1 

h > 



Ni ll c lli 

Taking expectations on both sides in combination with Prop. B.2. in [16] yields the desired 
result. ■ 
If L is large with uniform PDP the bound approaches the high SNR AWGN capacity and is 
thus too optimistic in general. 

C. Convergence to the ergodic capacity 

Let us now treat the ergodic case. The ergodic capacity is given by [4]: 

C e = E h (max {log , 0}) where £ such that 
P* = E h (max^-^O}) 
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Note that for low SNR the first order term of the DLC is not bounded with respect to L. 
Hence, the ergodic capacity has the same property which is in accordance with results in [9]. A 
similar convergence can be shown for the high SNR case. 

Corollary 4: Under the assumptions of Lemma 3 the DLC converges to the ergodic capacity 
as K — > oo. 

Proof: We only have to show that C e scales as log (P*) + H F as P* — > oo. We can again 
use a truncation argument. Let h\ := max {hi, e}. Then, we have for sufficiently large P* 



£-E h ((hl)- 1 )=P* 



and 



C e = E h (log(M) 

= E h (log((E h ((^ 1 )- 1 )+P*) h{)) 

= log (P*) + E h (log (h\)) + log (E h ((hi)' 1 ) (PT 1 + 1) 

Now, again by a bounded convergence argument 



limE h (log (hi)) =E h (log (hi)) 
40 



and 



C e 



log (P*) + H F 



1, P* -> oo, 



which completes the proof. ■ 
We have the following coding implications: 

V. OFDM BROADCAST CHANNEL DLC REGION 
The delay-limited region is defined as follows: 

Definition 2: A rate vector R* lies in the DLC region Cdl (P*) with sum power constraints 
P* constraints if and only if for any fading state h there is P' solving 

minP s.t. R*eC BC (h,P) 

and 

E h (P') < P* 
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Furthermore, R* is on the boundary £> DL (P*) if and only if 

E h (P') = P*. 

To evaluate the OFDM delay-limited region turns out to be very difficult. This is due to the 
fact, that we have only an implicit characterization for the DLC. This means, that we can check 
for every rate vector R, whether it lies inside the DLC region or not, simply by solving the dual 
problem for each h. To evaluate the expectation necessary to check the average power condition 
can be done by Monte-Carlo runs. However, it is very difficult to determine all rate vectors, 
which can be achieved with a fixed sum power constraint P*. 

Nevertheless and although computationally demanding, the OFDM-DLC region can be 
calculated up to any desired finite accuracy. To this end we restate the Algorithm 1 from [10] 
and define 



e n> m \2 + J2vh 12 {e nkU) ' k ~ 1) e»= i+1 \ (39) 

l \' l 7T k (m),k\ . =1 \ n "K k {j),k\ J J 

= [log(/i m ) + ™m,fc] • (40) 

Then Algorithm 1 yields the minimum sum power necessary to support a set of rates R. 

Algorithm 1 Iterative "Rate Water-Filling" 

Set R m)k = Vm e Ai, k = l, ...,K _ 

while desired accuracy is not reached do 
for m = 1 to M do 

(1) Compute the coefficients n m ^ (39) for user m 

(2) Do water-filling with respect to the rates R m ,k for user m according to equation (40) 
end for 

end while 



To evaluate the region Cdl (P*), first the single user DLC rates have to be calculated for all 
m. This is done by the evaluation of (7) for fixed Cdl (P*) and bisection, since Cdl (P*) is 
monotone in P*. Due to the convexity of Cdl (P*), any convex combination R int must lie inside 
Cdl (P*)- On the other hand, the single user rates R m are a component- wise upper bound for 
all other rate vectors. Since the necessary power P (aR^) is monotone in a, simple bisection 
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can determine the boundary of the region for each angle. For any new points on the boundary, 
the refinement procedure can be repeated until the desired number of points defining the border 
is obtained. This procedure is summarized in Algorithm 2. Note, that if the distribution and 
number of taps is the same for all users, the region is symmetric and can be constructed by 
mirroring one calculated sector. 

Algorithm 2 OFDM DLC Region Algorithm 

(1) Determine single user DLC (axis points) by evaluation of single user DLC (see below) 

and bisection 

while desired accuracy not reached do 

(2) For any two neighboring vectors Ri G £>dl (P*) and R 2 G £>dl (P*) on the boundary 
calculate interpolated vector H int = 1/2 (R x + R 2 ) 

(3) Adjust a > 1 by bisection using Alg. 1 such that «R illt G £>dl (P*) 
end while 



VI. OFDMA ACHIEVABLE DELAY LIMITED RATE REGION 

Compared to the single user case the multiuser DLC region is more difficult to analyze. In 
order to get some insight in this case we derive simple resource allocation schemes based on 
OFDMA and rate water- filling. To this end, we assume independence of the subcarriers achieved 
by complex Gaussian distributed path gains and uniform PDP. It is worth noting that this can 
be imagined as that we take only L independent frequency samples and assume that the other 
value are approximately equal in the neighborhood of the L subcarriers. Although seeming quite 
complicated, the following lemma yields a simple lower bound on the OFDM DLC region, 
implying an OFDMA strategy. 

Lemma 4: Let i = (h ■ ■ ■ , ,ik) G [1, M] K be a multi-index and let the set K} m count the 
number of user indices in i equal to m. Let X s C [1,M] K be the subset that contains these 
multi-indices, where all users occur in the multi-index i at least s times, i.e. 

X s = {\:K} m >s,meM} . 

Then the average required power P* to support any rate vector R in each fading state is upper 
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bounded by 



M 



m I R 1 

g m,k[p] — \ 



i S X s m=l p=l Sm -P 



M |^/ M J ( e R m,k[p] — 1 

^ + (M* -|X,I)E E 

m=l p=l 



M K e~ 1 

v m,p 



Where 



fc[p]e^ 
is the solution to 



c _1 



A 



AT 1 



0; E ^"i,fc[p] — 



p=l 



and 



f°° 1 



(41) 



with h' kh] being the p-th order statistic of |/C^| random variables h' with cumulative density 

F h , {x) = (1 - e~*) M . 



l k[ P ] 
function 



R 



■m,k\p] 



is the solution to 



pRm,k\p\ 



and 







\K/M\ 
p=l 



x fc[p] 



with /i'^pj being the p-th order statistic of [||J random variables h" with cumulative density 
function 



F h >. (x) 



M 



M-2 



x \ —x' 



dx'. 



l-e~ x 

Proof: The basic idea is to distinguish between two cases: The case, where each user has 
the best channel at least on s subcarriers and the case where at least one user has on less than 
s of the K subcarriers the best channel. With the multi-index i = (iy . . . , , i K ) e [1, M\ K define 
the event 

Wi := [uj : hi ul (u) > h hl (u) h iK>1 > h hK {u)\ mK } • 
Note that by the absolute continuous fading distribution, we have 

£Pr(W,) = l 
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since the remaining events occur with zero probability. Thus, we can express the average power 
P* by 

where I s C [1, M]^ is the subset that contains the elements where all users occur in the multi- 
index i at least s times. Let !O m be the set that counts the number of user indices in i equal to 
m. Fixing i and m and ordering the values according to 

the first term on the right hand side is bounded by 



^E h (p'i^)<EE E E * 

iel a m=l p=l 



n, (42) 



with -R^ fe [ p ] such that the required rates are supported: Since the expectation on the RHS is 
independent of the actual referred subindex in the multi-index i and depends only on the number 
of emerging entries of user m in i counted by the set K? m we can replace the RHS and rewrite 
the inequality in (42) as 



M K n \ h 

iSX s i£l s m=l p=l 



m,fc[p] 



e m - fc w — 1 

iel 3 m=l p=l Sm >P 

where ( m ^ p is the expectation of the p-th inverse channel coefficient. The distribution of the p-th 
order of the best channels /i™ ax for some k is independent of k and given in equation (9). The 
cdf and pdf / (x) of /i™ ax in turn can also be derived by the order statistic from (9) and be 
expressed as 

Pr > x, h 1A > h w ^) 



(x) 



Pr (h lt i > h w ^i) 
= M U-e - *) e _s tte 7 . (43) 

J X 



and 



/ i?) = = M (l — e ) e x (44) 
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Substituting (43) and (44) into (9) yields the desired densities. 

For the second term we need to carry out a different strategy. For all cases represented in the 
complementary set X s at least one user occurs in the multi-index i less than s times and hence 
has the best channel on less than s subcarriers. For the case s = 1, the strategy of the previous 
term can not even guarantee his delay limited rate requirement. Alternatively, we simply divide 
the set of subcarriers in M sub-bands where to each user [K/M\ subcarriers are allocated and 
do rate water-filling as done for the first term and take the best out of this set. Hence the second 
term is upper bounded by 



/ M [K/M\ 5 

£ Eh( pw)<5> E E 1 T^7 i 

M [K/M\ 



where the set HY 1 is defined as 



ivi L IS -/ lvl i / 

Yl { eRmMp] - l ) Eh ( 

rr> — 1 n— 1 ^ 



h 



m,k[p] 



Ti™ (45) 



n ? ■= > WHL^.-'^aH > ^aHL^[ ( 46 ) 

Note, that the second inequality stems from the fact, that the expectation is conditioned on the 
set H™, assuming that user m has not the best channel on any of his subcarriers. 

Since all subcarriers are independent, we define for each subcarrier the following conditioned 
probability and get after some manipulations 



F h „ (x) = 1 - Pr ( h rn<k > x\ h n<k > h rn<k , n ^ m) 

_ 1 _ Pr > x,h 2 ,i > Kim) 
Pr (h 2 ,\ > h w -t 2 ) 

= 1 - M J. ( 1 - e ") 

Thus, with (9) we can express the condition expectation as 



(47) 



Ey 



hm,k\p] 

leading to (41). Since the addends do not depend on the index i, the first sum in (45) can be 
substituted by the factor \X S | = M K — \T S \ leading to 

M \K/M\ f e R m . k [p] — 1 ) 

E Eh (P'\ WO < (M K - 127-1) E E • 

ieJs m=l p=l m >P 
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This concludes the proof. ■ 

Note, that for the case M — 2, the expression for the cdf in (47) simplifies to F h n (x) = 
1 — e~ 2x . Instead of partitioning the subcarriers equally among the users, it is possible to share 
them in any other relation. Then the second sum of the second term in (4) is not has not [K/M\ 
addends but K m addends for each user with J2m=i K m = K. So it is especially reasonable to 
share the subcarriers proportional to the users rate requirements such that K m = R m / Ylm=i Km- 

The bounds from the previous section are illustrated in Fig. 8 and Fig. 7. The bounds are 
shown for different values of s, changing the relation between term 1 and term 2 in (4). Fig. 7 
depicts the low SNR case. It can be seen that the lower bound nearly achieves the entire region. 
This is due to the fact that in the low SNR regime only the best subcarrier is used. This is 
realized with rate water-filling, even if not perfect but only ordinal information, e.g. the ranking 
of the subcarriers, is present. The remaining gap stems from the second term and the fact, that 
users can "collide", i.e. have a common best subcarrier. 

In contrast, in Fig. 8 the high SNR scenario is presented. The bound improves as s is increased 
up to s = 4. From s = 5 on, the bound degrades once again. This is since it is not optimal 
to support the entire rate only on one subcarrier, even if a user has only one best subcarrier. 
Thus, the bound improves as s is increased. Note, that for the case that both users have similar 
rate requirements, i.e. the sum DLC case, the bound achieves a major part of the DLC and 
outperforms the time-sharing strategy. The remaining gap on the axes is much bigger. However, 
the gap on the axes can be reduced by sharing the subcarriers proportional to the users rate 
requirements and using only the second term (thus making the conditioning of the pdf needless). 
This is illustrated with the curve called prorated FDMA. The discontinuity stems from the 
switching of the subcarrier allocation, since this is a discrete procedure. The dashed blue curve 
indicates an achievable OFDMA DLC region, which is given in Algorithm 3: For any rate vector 
R and any channel realization h the sum power minimization algorithm is usd to calculate the 
optimal resource allocation. The resulting Lagrangian multipliers /x are taken to allocate the 
subcarriers according to the maximum weighted channel rule m k = argmax me _vi ^wh m ,k- Once, 
the subcarrier allocation is done, the optimal resource allocation is obtained by water-filling such 
that the rate requirements are met. 

This scheme requires perfect CSI but is computationally still relatively simple due to the 
iterative water-filling principle. It can be seen from Fig. 8 that the algorithm yields good results 
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Algorithm 3 OFDMA DLC Algorithm 

(1) for given rate vector R and channel realization h solve the minimum sum power problem 

with Alg 1 

(2) assign subcarriers according to m k = argmax meA 4 ^ m h m ,k 
for m = 1 to M do 

(3) do water-filling with respect to the rates R m ^ for user m such that rate requirement R m 
is met 
end for 



where any other FDMA approach has to compete with. 

VII. Conclusions 

We studied the delay limited capacity of OFDM systems. We have shown that explicite 
expressions can be found in the low and high SNR regime even for the challenging correlation 
structure of OFDM. Even though we presented our results in the context of OFDM they are not 
restricted to this class but apply to other channels such as MIMO as well. On the other hand, 
still a basic open problem is the complete characterization of the DLC for all SNR and arbitrary 
power delay profile. Here, we were not able to give universal bounds and it is an interesting 
problem to show that the dependence is in general so-called Schur-concave which implies that 
a uniform profile maximizes the DLC in all cases. Furthermore, we analyzed the OFDM BC 
DLC region and derived lower bounds based on rate water-filling. In the low SNR regime and 
concerning DLC throughput, these bounds perform very well. To approach the DLC close to 
the axes in the high SNR regime, a prorated strategy has to be used. All bounds merely use 
order statistics and involve only ordinal - and thus partial - channel knowledge, which suggests 
savings for the design of future feedback protocols. Moreover, an additional FDMA strategy 
using full channel state information is proposed, performing very well over the entire region. 

Appendix 

A. Proof of Theorem 1 

By Theorem 3 in [13] we have to check that: 
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E h (exp 



e 4 



■+Tr i=z a i uj i +0{LJ 6 



holds for any 9 k 



2irk/K, k = 0, . . . , K — 1, and for any real number uo in the non-empty interval [—d, d] for 
some d > 0, and furthermore 



£ ES=i EJU*^ Eh (exp [#te (Ef =1 Q h^" 1 + 



= e-^-^+EU^dwil+I^D'+oCdwil+l^l) 6 ) holds for all real numbers uj u uj 2 in [-d,d] 2 . 

We show only the more complicated second condition. The first condition can be easily 
deduced by our assumptions on the distributions and observing that the subcarriers' real and 
imaginary parts are independent. We have 

L 



E h exp 



(W fcl )+u, 2 sin(W fc2 )) o((M+M) 3 ) 



By analytic expansion of the first factor and observing that 

K-l K-l L 

AL 



^ K-l K-l L 

7T^2 ^ (^i cos (W k J +uj 2 cos (W k2 )) 2 + sin (W kl ) +u 2 sin (W k2 )Y 



ki=0 fc 2 =0,fc 2 ^fci 1=1 



K-l K-l 



- + -j- + 2ujiuj 2 ^2 ^2 S ( cos (^ fe i ) cos (^ fe 2 ) + sin (^ fe i ) sin ) ) 



k 1 =0k 2 =0,k 2 =£k 1 1=1 



4 4 



where the last step is due to the standard trigonometric relation 



K-l 



cos (w k 



k=0 



K 1 = 
1^0,1<K 



we have finally 



Tl E; L =i(^icos(W fel )+aj 2 cos(«0 fc2 )) 2 + (a;i sm(i6» fcl )+w 2 sin(Wfc 2 )) 2 +0((|wi |+|w 2 1) 3 ) ) 



J l "2 



e 4 4 



+£L3^(M+Mr+o((M+M) 6 ) 



which is the desired result. The proof that non-uniform PDP can not improve this bound follows 
from union bound and is omitted [13]. 
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B. Proof of Lemma 2 

The order of the lower and upper bound was already derived in [13]. Due to the correlation 
structure imposed by the FFT and Rayleigh fading with uniform PDP the channel gains 
fyfci+(o-i)fc 2 > ^2 £ {1, . . . , L}, are independent for any k\ G {1, . . . , a} for some a e N where 
a = K/L. Since the maximum || (h kl , h kl+a , . . . , h kl+a (L-i)) 1 1 is below or equal x if /loo is 
below or equal x but not conversely, it follows the lower bound [13] 

Pr < x) < e~ L£ 

which is the desired lower bound if we set e = 7 log [log (L)] / log (L) (in fact, obviously, it is 
even stronger and can be to strengthened to fall off with L^ 1 instead of order log -1 (L)). The 
upper bound is obtained by observing that for any a > 1 we have by the FFT structure 

PrO^ > x) < a \l- (l-e~ x y 

Setting this time x = (1 + e) log (L) for some e > yields 



Pr(/ ioo >(! + £) log (L)) < a 



1-1- 



all — exp 



< a ^1 — exp 

aL~ £ 

< 
< 



1 _ L-a+ £ ) 



1 _ 



using log (1 — x) > —j^; and e 31 > 1 — x. Hence, if we set e = 7 '"log"! we nave 



Pr (/ioo > log (L) + log [log(L)])< 



K 



log 7 (L) ' 

Combining this with the stronger lower bound yields the result. 



C. Proof of Proposition 7 

We can frequently use the following well-known inequality: suppose that fi, . . . , f K are 
functions defined on a domain f2 equipped with some measure F with f k e £ Pfc (fi) and 
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E^P^lthen [17]: 



[f[fk (x) dF(x) <f[( f f Pk (x) dF (x)) 



(48) 



This inequality is tailored for the situation at hand: suppose that for some a,b > 





b 



dF h (h) 



which yields by application of (48): 



a 



fc 1= l l JK + k 2 =l 



'fcl + (fc2 — l)l 



dF h (h) 



(49) 



The inner term on the RHS of (49) generally means multidimensional integration with usually 
dependent random variables which cannot be directly carried out. Hence, we resort to some 
bounding techniques and have to show that an upper bound holds for the inner integral 
independent of k±. 

In order to obtain an upper bound on the inner term on the RHS of (49) choose some e > 
and subdivide the integration domain in parts where in each dimension the range of integration is 
either in the interval [0, e] or outside this interval. For those dimensions that are within this interval 
we bound the corresponding marginal density by the constant c s and calculate the remaining 
integral while for those dimensions that are outside the interval we simply set the values of the 
integrand to e. Suppose that I dimensions are within the interval. Since it does not matter what 
particular dimensions are chosen for this decomposition we have I out of b possibilities that can 
be equally treated. Hence, we can write 




b 



-+ k=l 
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and since 

E ©a -«>'«" = ! 

for any < e < 1 it follows 

E h (h) < Cs b(^-^j +(l-c a )b 

where e = 1/6 < 1. 

D. Proof of Proposition 8 

We can nicely use Holder's inequality. Swapping expectation and product operator we have 
by partial integration 

/ h^f{h)dh=h- b +1 f{h) -—r h- b +1 f{h)dh 



=0 

6 



-/ h~ l b +1 f (h) dh 
1 7r + 



since / finite everywhere. Define := {h : /' (h) < 0} and let < a i <b i < 00, i G N, be 
the interval boundaries of /i where /' (/i) < as well as < a\ < bf < 00, % e N, those where 
/'(/i) > 0. Representing (-/') as (-/') = (-/')^ ° (-/')* in R + yields 



- / h~i +1 f (h) dh=-±- I h*? (-/' (fc))^ (-/' (/*))* dh 

1 JRI 0—1 Jjjo 



b 

Setting p — b/ (b — 1) and g = 6 



/ h-* (-f(h))dh<[ f (-f(h))dh) If h(-f(h))dh 



Again by partial integration for the last term 



hf'(h)dh = -hf(h)\™+ [ f(h)dh 



< 1 



38 



and since 



— T / „ h~i+' f (A) dh < -£ f f (h) dh 

^12"? (f (K) - f (4)) 



" ■. ... 



we have finally 

/ h-*f(h)dh 



< 



b 



1/6 b 



i>l / i>l 



6-1 

which concludes the proof. 



E. Proof of Proposition 9 

Suppose K, L to be even integers. We can apply inequality (49) with a = K/2 and 6 = 2: 

n dp * ( h ) ^ n / K h i i/2k w (h) 

+ fe=i i=i V/ M £ / 

Since h k = \c T C k \ 2 with C k ■= [l,e-^ k - i y K ,...,e~ 2ir ^ k - 1 ^ K - 1 ^ K ] T we can write for the 
inner integral 

+ J J R2L k=l k=l 

where f r , f l are the bounded densities of real and imaginary parts of the complex path gains and 
(•) r ' 1 is a shorthand notation for real and imaginary part operators. Let I e [1, K/2] be arbitrary 
but fixed. The following change of coordinates is based on the observation that for K, L even 
Ch Ci+k/2 are orthogonal in C L , VZ. For / = 1 this is obvious since the first and K/2-th vector 
consist of an even number of binary ±l's only. For I > 1 the same follows from the fact that 
two orthogonal vectors remain orthogonal if the are both multiplied by the same complex phase 
factors. 

Hence, there exist Ct , i — 2, • • • L — 2, that can be chosen to span the basis of the orthogonal 
complement. By change of coordinates c — > ^h 2 j ~h\ = c T C, h h z 2 = c T C K / 2+ i,h z = c T C, z i) i = 
2, ... L — 2, where Ci+k/2: Ct > is an orthogonal transformation and the Jacobian equals 
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1/L L and, further, h% — > (\h z k \ 2 , ip z k ) = (R z k , ip%) , k — 1, . . . L, of which the Jacobian is 1/2 \/k 
we obtain by assuming /gr (x) < i>^/ 2 e~ a l x l 2 and /gi (x) < i;^ 2 e -a l x l 2 , Vfc: 



RHS of (50) < ^ / 



fc=i 



-i—r f i afl l T—r /" aR k 



TL 

rL / r \ L-l 



< rT - [ — ) UJU 1 P^. 8) 



Remark 4: For many fading distributions the claim f-r,i (x) < v l J 2 e~ a ^ 2 might be too 

k 

restrictive and shall be replaced by f 5 r,i(x) < max |c (max) , v]J 2 e~ a ^ 2 1 where c (max) > is 
some global constant. The latter inequality, however, does not separate over R L as required in 
the proof here. In this situation, we can obtain better bounds for specific cases by combining 
the techniques of Prop. 7 - Prop. 9, e.g. by splitting up the integration domain similar to Prop. 
7. 
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Fig. 1. OFDM DLC and lower bounds for L=K=16. 




Fig. 2. Scaling in the low SNR region: The dashed lines indicate the scaling at low SNR given by So = l/A"log 2 (l + 
KP* \og{K)) 
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Fig. 3. Scaling in the high SNR region: The black line indicates the scaling at high SNR given by S = log (P*) + Hf (for 
Rayleigh fading). The dashed lines give the OFDM DLC for L = K = 2, 4, 8, 16. 




Fig. 4. Comparison of the DLC with the first and second order approximations and the lower bound from Fig. 2. Channel 
with 64 taps delay spread and system with 64 subcarriers 
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Fig. 5. Comparison of the DLC with the first and second order approximations and the lower bound from Fig. 2. Channel 
with 1024 taps delay spread and system with 1024 subcarriers 




0.5 1 1.5 2 2.5 3 3.5 



R 1 [bps/Hz] 



Fig. 6. Example for an iteration of the described algorithm to calculate the OFDM DLC region. OFDM MAC DLC region for 
2 users with 7 i.i.d taps each and 16 subcarriers at WdB 
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Fig. 7. OFDM DLC region for 16 subcarriers at -20 dB. The lower bound is shown for different values of the parameter s. 
The lower bound degrades for increasing values of s. 




Fig. 8. OFDM DLC region for 16 subcarriers at 10 dB. The lower bound is shown for values of the parameter s = 1, ...,4. 
Further, the prorated FDMA region and an achievable OFDMA DLC region are depicted. 



